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An  interesting  electromagnetic  radiation  mechanism  known  as  the 
electron  cyclotron  maser  has  been  the  subject  of  intense  research  activities 
/ in  recent  years.  This  mechanism  has  been  the  basis  for  a new  class  of 

microwave  devices  called  gyrotrons  capable  of  generating  microwaves  at 
unprecedented  power  levels  at  millimeter  and  submillimeter  wavelengths. 

A detailed  description  of  the  cyclotron  maser  mechanism  is  given  in 
Reference  U,  and  brief  summaries  of  gyrotron  theories  and  experiments 

5.6 

together  with  lists  of  references  can  be  found  in  recent  review  papers. 

In  the  present  study,  we  will  concentrate  on  a particular  type  of 

gyrotrons— the  gyrotron  travelling  wave  amplifier  (gyro- TWA).  The  basic 

physical  processes  taking  place  in  a gyro- TWA  have  been  analyzed  in  recent 

U 7-12 

linear  and  nonlinear  theories.  * In  the  actual  operation  of  a gyro- TWA, 

the  beam-to-wave  energy  conversion  efficiency  Is  perhaps  the  most  important 
consideration.  References  10-13  contain  detailed  studies  of  the  saturation 

« 

mechanisms  and  calculations  of  efficiency  for  the  operation  at  the  fundamental 
cyclotron  harmonic.  However,  the  scaling  of  the  efficiency  with  respect 
to  the  various  modes  and  parameters  of  operation  has  not  been  considered 
in  any  detail,  nor  has  the  operation  at  the  nonfundamental  cyclotron 
harmonics.  The  nonfundamental  harmonic  operation  is  of. great  importance 
for  the  amplification  of  submillimeter  waves  if  unrealistically  high  magnetic 
fields  (>  100  kG)  are  to  be  avoided.  In  anticipation  of  the  growing 
experimental  effort  aimed  at  the  generation  of  submillimeter  waves,  our 
main  purpose  in  the  present  study  is  to  derive  a general  analytical  expression 
* for  the  operating  efficiency. 

Derivation  of  a general  dispersion  relation.  The  typical  configuration 
of  a gyro- TWA  consists  of  an  annular  electron  beam  propagating  inside  a 
waveguide  of  circular  cross-section  (of  radius  rw).  The  electrons,  guided 

Note:  Manuscript  suomitted  March  22,  1979. 
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by  a uniform  magnetic  field  (Bq  g^),  move  along  helical  trajectories. 

In  our  model,  all  the  beam  electrons  have  the  same  perpendicular  velocity 

(v  ),  and  also  the  same  parallel  velocity  (v  ),  with  their  guiding  centers  t 

jlO  zo 

uniformly  distributed  on  a surface  of  constant  radius  rQ  (see  Figure  1). 

We  assume  that  the  beam  is  sufficiently  tenuous  that  its  space  charge 

electric  field  can  be  neglected,  and  the  spatial  structure  of  the  vacuum 

waveguide  mode  is  unaffected  by  the  presence  of  the  beam.  The  beam 

interacts  with  a single  TE^  waveguide  mode,  where  m and  n are,  respectively, 

the  azimuthal  and  radial  eigenmode  numbers.  The  dynamics  of  the  electron 

beam  is  described  by  the  linearized  relativistic  Vlasov  equation  which 

couples  with  the  Maxwell  equations  to  form  a complete  set.  Methods  for 

obtaining  the  dispersion  relation  from  this  set  of  coupled  equations  are 
7 

standard.  Here  we  present  the  result  directly. 
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where  0^  = v Jc,  0C  = eB^mc,  rL  = v is  the  n-th  root  of 

J^(x)  = 0,  Jm(x)  is  the  Bessel  function  of  order  m,  J^(x)  = dJ^(x)/dx, 


2 2 

v s Ne  /me  is  a dimensionless  beam  density  parameter,  N is  the  number  of 
electrons  per  unit  axial  length,  and  the  functions  k , H , 0 are 

mr> 9 SHi  iSm 

defined  as  follows:  K_  = ^(x )-J  , (x  )j  . (x  ),  H (x.y)  = 

tan  mv  mn'  m-1'  mn  m+lv  mn  * smv  ,jr/ 
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[Js-m(x)J3(yi]2’  and  Qsm(x>y)  = 2Hsm(x’y)  + y [JLn(x)Js(y)  Js(y)  + l^-m-l(x) 
Jg(y)Jg_l(y)  - 11,6  disPersion  relation 

In  Equation  (1)  is  more  general  than  those  found  In  literature  In  that  the 
waveguide  mode  numbers  m,  n,  and  the  cyclotron  harmonic  number  s are  all 


treated  as  free  parameters. 


This  will  provide  the  option  to  analyze 


all  possible  Interactions  as  one  searches  for  the  optimum  efficiency.  On 

the  right-hand  side  of  Equation  (1) , the  first  term  is  the  source  of  the 

instability,  while  the  second  term  imposes  a threshold  beam  energy  for  the 

instability.  For  the  non fundamental  harmonics  (s>l),  it  can  be  shown  that  the 

threshold  energy  (typically  below  l keV)  is  much  lower  than  the  typical 

beam  energy  (tens  of  keV);  hence  the  second  term  can  be  neglected. 

Substituting  cu  = U)  + Au>,  k = k into  Equation  (1),  where  (u>  ,k  ) 

o z zo  o zo 

is  the  point  at  which  the  waveguide  characteristic  curve 

2 ,2  2 22.2 

u>  - k c - x c /r  *=  0 (2) 

z mn  v 

intersects  with  the  beam  characteristic  curve 

<"  - Vzo'  = 0 (3) 

we  can  easily  evaluate  Atu(»Au)r  + lAu^)  , with  the  result 
r.  2 „ 2 411/3 


2 2 1 
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mn  smBlo 
— 
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In  Eqatlon  (4),  Aa)^  gives  the  frequency  width  for  resonant  beam-wave  inter- 
action, and  Atu^  gives  the  growth  rate. 

Numerical  simulations.  The  preceding  linear  analysis  has  been 

13 

complemented  by  a single-wave  numerical  simulation  code  developed  to 
simulate  the  amplification  of  azimuthally  symmetric  modes  at  the  nonfundamental 
harmonic  frequency.  The  simulation  shows  two  different  but  simultaneously 
present  saturation  mechanisms- -depletion  of  free  electron  energy  and  loss 


r-J-  . 


of  phase  synchronism.  The  first  mechanism  dominates  when  the  beam  energy  is 
slightly  above  threshold.  Saturation  occurs  as  soon  as  the  beam  loses  a small 
amount  of  energy  and  the  system  becomes  linearly  stable.  The  second  mechanism 
dominates  when  the  beam  energy  is  well  above  the  threshold.  Saturation  occurs 
because  an  average  electron  loses  so  much  energy  that  its  relativistic 
cyclotron  frequency  no  longer  matches  the  wave  frequency  to  favor  unstable 
interactions.  Figure  2 shows  a typical  efficiency  curve  as  a function  of 
time.  It  increases  at  the  linear  growth  rate,  then  exhibits  an  oscillatory 
behavior  after  saturation.  As  founa  in  Reference  11,  both  mechanisms  are 
important  for  the  fundamental  cyclotron  harmonic  interaction.  However,  there 
is  a basic  difference  between  the  fundamental  and  nonfundamental  harmonic 
interactions;  namely,  for  the  latter  interaction,  we  find  that  the  threshold 
beam  energy  is  typically  so  low  that  the  energy  depletion  saturation  mechanism 
can  be  disregarded.  This  has  allowed  us  to  derive  an  analytical  scaling 
relation  for  the  efficiency. 

Efficiency  scaling  relation.  For  the  purpose  of  efficiency  optimiza- 
tion in  the  multiple  parameter  space  typical  of  a gyro-TWA>system  an  analytical 
scaling  relation  is  essential.  Without  loss  of  generality,  we  shall  derive  the 
efficiency  scaling  relation  in  the  beam  reference  frame  (i.e.,  the  frame  in 

which  vzo  = 0) , and  denote  the  beam  frame  quantities  with  a prime.  A simple 

11  12 

Lorentz  transformation  ’ can  be  used  to  convert  the  beam  frame  efficiency 


into  lab  frame  efficiency.  In  order  to  minimize  the  possibility  of  spurious 
oscillations,  it  is  advantageous  to  choose  a magnetic  field  such  that  Eqs.  (2) 
and  (3)  Intersect  at  only  one  point  (i.e.  at  a grazing  angle).  In  the  beam 


frame,  this  Implies 


<t  s s n = x c/r 
O c mn  w 


-4- 


Henceforth,  our  analysis  will  be  restricted  to  this  particular  case. 


In  the  linear  analysis,  we  have  shown  that  the  following  condition 
holds  at  the  onset  of  the  instability, 

I 

A u) 


u,  . - 8 0 /v 

r e o 


(6) 


where  Qg  ■ eBo/mc  and  Au>r  is  a positive  quantity  given  by  Equation  (4).  As 
the  instability  evolves,  the  average  energy  of  the  electrons  decreases,  hence 
the  left-hand  side  of  Equation  (6)  also  decreases.  Knowing  this  tendency 
and  the  width  for  resonant  interaction,  we  expect  the  following  condition  to 
hold  at  the  saturation  state, 

(7) 


(!) 


s 0 /<y  > -Au/ 
r e s r 

i I 

where  <YS>  is  the  average  y of  all  the  electrons  at  saturation.  This  is 
the  physically  expected  saturation  condition  because  any  further  decrease  in 
electron  energy  would  shift  the  electron  cyclotron  frequency  out  of  the  range 
for  resonant  beam-wave  interaction.  Equation  (7)  was  found  to  be  in  good 
agreement  with  our  extensive  simulation  data. 

From  Equations  (4)  - (7)  and  assuming  A* '«s)^ , we  obtain  an  approximate 

• t2. 


estimate  for  the  efficiency  T|, 

• Yo_<Ys>  i 
7)  ==  -°r  ~ ~ 

Y -1 


V1 


1* 


2 v*Yo  H (x  r /r  ,x  r.  /r  ) 
sm  pm  o v inn  L w 


K x 
mn  inn 


(8) 


where  H , K , and  x are  all  frame  independent  quantities, 
sm*  mn’  mn 

Equation  (8)  shows  how  T)  scales  with  the  mode  numbers  m,  n,  s,  and 

it  i 

the  operating  parameters,  , yQ  , and  rQ.  Note  that  H is  independent  of  the 

waveguide  radius  v^.  The  scaling  of  7) * with  respect  to  v'  is  especially  simple 

namely,  For  a given  mode  of  operation,  the  choice  of  r^  should  be 

such  that  H falls  on  or  near  its  peak  value.  Scaling  of  7} ' with  respect  to 
sm 

m,  n,  s,  and  V0  is  more  complicated  and  can  best  be  seen  through  numerical 
plots  of  Equation  (8).  Some  examples  are  shown  in  Figure  3,  in  which  TJ  is 


is  plotted  (in  continuous  curves)  against  the  beam  energy  (y©)  for  m“0, 

3=01*2,  3,  and  4.  In  each  case  the  beam  position  (rQ)  has  been  chosen  to  maximize 
Han.  It  is  seen  that  higher  cyclotron  harmonics  generally  give  lower  efficiency, 

i 

as  expected,  and  the  behavior  of  H as  a function  of  the  beam  energy  is 
different  for  different  harmonics.  For  comparison,  the  numerical  simulation 
data  was  also  shown  on  the  same  figure  (dots).  Since  11  is  an  oscillating 
function  of  time  (see  Figure  2),  we  have  consistently  chosen  from  the  simula- 
tion data  a 11'  which  is  80  percent  of  the  peak  value.  The  agreement  between 
Equation  (8)  and  the  simulation  data  is  very  good  for  all  the  non fundamental 
harmonics  except  the  second  where  there  is  a small  discrepancy  at  lower  beam 
energies.  This  is  due  to  the  fact  that  the  energy  depletion  saturation 
mechanism,  which  still  plays  a minor  role  at  the  second  harmonic,  has  been 
neglected  in  the  analytical  theory. 

To  summarize,  we  have  presented  a general  linear  dispersion  relation 
and  a concise  efficiency  scaling  relation  for  the  gyro-TWA.  These  results 
represent  the  first  comprehensive  analytical  and  simulation  studies  of  the 
non fundamental  harmonic  interactions.  In  terms  of  experimental  applications. 
Equation  (8)  can  be  used  to  asses  the  feasibility  of  a given  experimental  goal 
and  to  select  the  optimum  mode  and  operating  parameters  for  achieving  it,  while 
Equation  (1)  can  be  used  to  calculate  the  properties  of  a gyrot-TWA  system, 
such  as  the  bandwidth  and  amplification  rate,  etc. 
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Fig.  1 - Cross-sectional  view  of  the  gyro-TWA  model.  The 
applied  magnetic  field  (not  shown)  points  toward  the  reader. 
The  electrons  are  monoenergetic  and  all  have  the  same  Larmor 
radius  r . Guiding  centers  of  all  electrons  are  uniformly 
distributed  on  the  circle  of  constant  radius  r . 


